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LANCZOS ALGORITHM - FE VERSION 
rj = K-'M qj 
r j =7 i - “jOj - Pft-i 

“j = q/Mfj 

Pj = (r j -i T Mr J _,) l/2 

Qj+1 = d/p j+ i)rj 
P j+1 = (r j T Mr j ) 1 ' 2 


The algorithm for recursively forming Lanczos vector qj tl from Lanczos vector q 

for the case of nonsingular K is shown. It involves an orthogonalization step and a nor- 
malization step. 


STARTING VECTOR 
F = be(t) 
r 0 = K" 1 b 
qi = (l/pi)r 0 
Pi = (r 0 T Mr 0 ) I/2 


Lanczos modes are most useful when the spatial distribution of the excitation is 
constant. This load distribution then determines the starting vector. 
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LANCZOS EQUATION FORMAT 


Mii + Ku = F 

MK -1 Mii + Mu = MK -1 F 

u = Qm z 

Qm = tqiq2— q m l 

Q m T MK-'MQ m i + Q m T MQ m z = QJmR-'F 

Lanczos vectors are used in a mode-superposition manner which is not exactly the 
same as the familiar Rayleigh-Ritz version of mode-superposition. Multiplication by the 
equation of MK "smooths" the loading. 

LANCZOS EQUATION FORMAT - CONT. 

Tm z + z = 8m 

g m = QJmk^f 

<*1 P2 
h a 2 P3 
P3 • • 

T m = 

• • Pm 

Pm«m 

The final Lanczos equations have a tridiagonal generalized mass matrix and have a 
unit matrix for the generalized stiffness matrix. The form of the generalized force vector 
g m is very special. 
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EXAMPLES 


• Lanczos Modes and Equations 

• Comparison of Normal Mode Models and 
Lanczos Mode Models 

+ Poles and Zeros 

+ Frequency Response Functions 

+ Transient Response 


U 1 U 2 U 3 u < u 5 u 6 

-4 - 4 - - 4 - 4 >- - 4 - - 4 - 

f(0 — | | | | 
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FINITE ELEMENT MODEL 


F(r) 


U ({) 



CONTINUOUS MODEL 


m\\ 


FE LANCZOS MODE MODEL 


T m z + z = g m 
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Note that there is tridiagonal inertia coupling of the Lanczos generalized coordi- 
nates, but note that the only Lanczos coordinate directly excited by the external force is 
the first coordinate. 
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SUN - MQ31 1 1 1 
1ST NORMAL VECTOR 


RUN - MQ3 1 1 1 1 
2ND NORMAL VECTOR 





2.00 3.00 

LOCATION 


RUN - M0311 1 1 
1ST LANC20S VECTOR 



RUN - W031111 
2ND LANCZOS VECTOR 



2.00 3.00 

LOCATION 


Typical normal modes and Lanczos vectors are shown here for axial deformation 
of a clamped-free rod. The starting Lanczos vector is based on a single force applied at 
the "free" end. Although finite element normal modes do have some stain near the free 
end, "exact" normal modes would all be strain-free at the end where the excitation force 
is applied. 




TRANSFER FUNCTIONS IN MODAL COORDINATES 


s 2 +2.48152 
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s 2 +70.87551 
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1.62980 
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TRANSFER FUNCTIONS - NORMAL MODE MODELS 
u,(s)/f(s) = 10.3415 1 [s 2 +3.28270) 2 ][s 2 +(6.82929) 2 ] 

[s 2 +(10.94259) 2 ]/A 4 

A 4 (s) = [s 2 +(l .57528) 2 ][s 2 +(4.83424) 2 ] 

[s 2 +(8.4 1 876) 2 ][s 2 +( 12.49644) 2 ] 

u,(s) /f(s) = 4.23934[s 2 +(3.52834) 2 ]/A 2 (s) 

A 2 (s) = [s 2 +( 1 .57528) 2 ] [s 2 +(4. 83424) 2 ] 

fl,(s) /f(s) = 2.02296 / [s 2 +(1.57528) 2 ] 


TRANSFER FUNCTIONS - LANCZOS MODE MODELS 
u,(s)/f(s) = 1 8.29893[s 2 +(3. 1 7763) 2 ][s 2 +(6.590 1 4) 2 J 

[s 2 +(10.97092) 2 ]/A 4 (s) 

A 4 (s) = [s 2 +( l .575285) 2 ][s 2 +(4.83427) 2 ] 

[s 2 +(8.48668) 2 ][s 2 +( 15.21 009) 2 ] 

u,(s) /f(s) = 7.14344[s 2 +(3.061 12) 2 ]/A 2 (s) 

A 2 (s) = [s 2 +( 1 .57529) 2 ][s 2 +(5. 1941 4) 2 ] 


u,(s) /f(s) = 2.51446 / [s 2 +(1.5857) 2 ] 


POLES AND ZEROS OF REDUCED-ORDER MODELS 


EXACT 

POLES 

(6DOF) 

4DOF 

NORMAL 

4DOF 

LANCZOS 

2DOF 

NORMAL 

2DOF 

LANCZOS 

1.57528 

1.57528 

1.57529 

1.57528 

1.57529 

4.83424 

4.83424 

4.83427 

4.83424 

5.19414 

8.41876 

8.41876 

8.48668 

- 

- 

12.49644 

12.49644 

15.21009 

- 

- 

TABLE 2. POLES OF REDUCED-ORDER MODELS 

EXACT 

POLES 

(6DOF) 

4DOF 

NORMAL 

4DOF 

LANCZOS 

2DOF 

NORMAL 

2DOF 

LANCZOS 

3.17759 

3.28270 

3.17763 

3.52834 

3.06112 

6.57266 

6.82929 

6.59014 

- 

- 

10.39230 

14.63693 

18.85315 

10.94259 

10.97092 

- 

- 


TABLE 3. ZEROS OF REDUCED-ORDER MODELS 


STATIC RESPONSE COMPARISON 


Normal 

Lanczos 

No. of 

Mode 

Mode 

Modes 

Model 

Model 


(fi/f)s=o 0-96961 

0.99951 

4 

" 0.91005 

0.99980 

2 

" 0.81521 

1.00000 

1 

' " 1.00000 

1.00000 

6(EXACT) 


TABLE 1. TRANSFER FUNCTIONS EVALUATED AT s=0 


These two tables compare poles and zeros of normal mode models and Lanczos 
vector models and the static response of normal mode and Lanczos vector models. 
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5. 00E-02 FREQUENCY (HZ) 5. OOE+UO 


A2: FRF (2 NORMAL) + + + 

A3i FRF (6 NORMAL) 

A4, FRF (4 NORMAL) o o o 


FREQRESP-BODE 

1 Y+ 1 Y+ # 3 

092485-001408 
092685-152824 



Ah FRF 


A3, FRF 


A4, FRF 


(6 LANCZOS) 

(4 LANCZOS) o o o 
(2 LANCZOS) 4- + + 


FREQRESP-BODE 

1 Y+ 1Y- n 4 
092485-002412 
021786-215523 


Direct frequency response functions (tip displacement/tip force) are shown for nor- 
mal mode models and Lanczos vector models. The Lanczos models have improved 
low-frequency performance. 
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Control maneuvers, such as slewing, are transient response problems. The next 
examples compare transient response solutions of normal mode an Lanczos vector 
models when the structure is subjected to step and ramp excitation. Two force distribu- 
tions are illustrated. 
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o 

o 


RUN 


MS1 1 1 12 


RUN - MD11112 




RUN - LD11112 



g 


RUN - LS1 1 1 1 2 


o 

m 


« 


o 



{ — 


f 


— 



s 


0.00 


2.00 3.00 

LOCATION 


e.oo 


These are modal and Lanczos solutions for a step force applied at the tip. Both 
displacement and strain solutions are shown. Note that the Lanczos solutions converge a 
little better than the normal mode solutions do. 
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DISPLACEMENT DISPLACEMENT 




X<"- 

\ N 


2.00 <S ® 4 JOO' 

LOCATION 


* 00 * 5 . 06 *' 


2.00 3.00 

LOCATION 


RUN - LD 1 1323 


RUN - LS1 1323 



2.00 3.00 

LOCATION 



2.00 3.00 

LOCATION 


These are modal and Lanczos solutions for opposing forces applied at nodes 1 and 
3 with a ramp time history. Note that the modal solutions are very poor, while the Lanc- 
zos solutions show excellent displacement and strain convergence. 




CONCLUSIONS AND RECOMMENDATIONS 


• Lanczos mode models represent low- 
frequency forced response better than do nor- 
mal mode models. 

• Lanczos mode models can be developed for 
both continuous and finite element structural 
representations. 

• Lanczos mode models for systems with multi- 
ple inputs and/or rigid body modes should be 
developed. 

• Numerical stability of the Lanczos algorithm 
should be assessed. 

• Control system designs employing Lanczos 
mode models should be attempted. 


